Optical properties of an interacting large polaron gas 

V. Cataudella"^, G. De Filippis*, G. ladonisi^ 

'^INFM and Dipartimento di Scienze Fisiche, Universitd di Napoli 1-80125 Napoli, Italy 
*INFM and Dipartimento di Scienze Fisiche, Universitd di Salerno 1-84-081 Baronissi (Salerno), 

Italy 

(February 1, 2008) 
The normal state conductivity, (t{uj), of a system of interacting large 

polarons is calculated within the Random Phase approximation and some 

numerical results are presented. The behaviour of the optical absorption as 

a function of the charge carrier density and of the temperature is analyzed 

for different values of the electron-phonon coupling constant. It is shown that 

(t{ijo) exhibits features similar to those observed in the infrared spectra of the 

cuprates. 
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I. INTRODUCTION. 



The infrared absorption of large polarons is a 



well studied problem.lll However, a large 



part of the studies has been focused mainly on the very low density regime (non interacting 
polaron limit) which is suitable for not heavily doped polar semiconductors and ionic insu- 
lators. In this regime the most accurate approach has been discussed by Devreese et al.i 
starting from the expression of the impedance function derived by Feynmann.i More recently 
the interest for this problem has been renewed in connection with the infrared absorption 
of cuprates in the normal phaseflJi In this materials the interaction among polarons can 
be important. It is, therefore, of interest to consider the effect of the Coulomb long range 
interaction on the large polaron absorption. To our knowledge this effect has not been yet 
considered within the Devreese-Feynman approach to the polaron absorption. 

The aim of this paper is to calculate the conductivity (j[u!) of a large polaron gas within 
the random phase approximation (RPA) by using, for the lowest order polarization insertion, 
the propagator of the Feynman polaron model. 10 We find that it is possible to identify in 
cr(co') different contributions with features common to some structures of the infrared spectra 
of cuprates. In particular we observe a significant displacement of spectral weight from higher 
to lower frequencies when the polaron density increases. 



We consider a system made of band electrons in the effective mass approximation in- 
teracting with non dispersive LO phonons and repelling each other through the Coulomb 
potential screened by the background high frequency dielectric constant BoM- The electron- 
phonon (e-ph) interaction is assumed to be: 



II. THE MODEL. 
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is the Frohlich e-ph matrix elemenlJ§. 

In the above expressions c^a (ct^) and (at) indicate, respectively, the annihilation 
(creation) operators for electrons and phonons, V is the system's volume, a is the Frohlich 
e-ph coupling constant, Rp = {%/ {2mujQ))^^'^ is the polaron radius and c^o is the LO phonon 
frequency. 

The normal state conductivity of this system, cr(c(j), is related to the retarded form of 
the current-current correlation function T{{uj) by the Kubo formula0: 



where n is the density of charge carriers. In order to estimate the conductivity we follow 
the approach suggested by Feynmani and adopted with success by Devreese et al.i for the 
optical absorption of non interacting large polarons. Within this approach the real part of 
the conductivity is given by: 

Re [a{u)] = -G{u) IjJ]^ 

where G[lS) = ne'^/m and ^{00) = II^{uj)uj/G{uj) (for the definition of n°(ci;) see below). A 
justification of eq.(^ has been given by Feynmani, while a more general derivation can be 
obtained by making use of the Mori-Zwanzig0 projection operator technique, as shown by 
Devreese et al-S and MaldagueEl. It is worth to note that the sum rule for the real part of 
the conductivity: 

r dcuRe [aico)] = ^ (5) 
JO 2m 

is satisfied by using eq.(|D. 

To evaluate Il^{uj) we use the set of diagrams shown in Fig.l, i.e., we consider a 
R.P.A.-like approximation for polarons interacting each other through the Coulomb po- 
tential screened by eoo {v^)- Following Mahan0, we write: 
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In the eqs-d^) and fi = x,y or z, e{q, u) is the dielectric function, M^is the renormahzed e- 
ph matrix element, D{q, uj) is the phonon Green function and ^^(ci;) is the boson occupation 
number. In our model the dielectric function of the system is assumed to be (see Fig.l): 



e{q,u;) = l-v^P{q,u;) 



(8) 



where P{q,u), the lowest order polarization insertion, can be expressed in terms of the 
polaron spectral weight function A{q,uj): 

h J (27r) "'-00 ZTT J-oo ZTT \UJ + S — U + to J 

Here np{u)) is the fermion occupation number. For simplicity we shall use in the Eq.(^ the 
unperturbed phonon Green function D^{uj) = 2ujq/ (cj^ — cUq + i5) and we shall ignore the 
frequency dependence of the dielectric function in the expression of the renormalized e-ph 
matrix element assuming: 
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Until this point we have not specified the polaron spectral weight function that appears in 
the Eq.(^). It is well known that, from a theoretical point of view, the problem of an electron 
in a polar crystal which interacts with the longitudinal optical modes of lattice vibrations has 
not been solved exactly. However, it is universally recognized that of all published theories 
that of Feynmani, which uses a variational method based on path integrals, gives the best 
available results in the entire range of the coupling constant a. For this reason we choose 
the spectral weight function of the Feynman polaron mo delii: 
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where = - l) jv, Mt = v^jw^, ^ = § [nuiv) (^^(f) + 1)]^/^ and h are the Bessel 
functions of complex argument. The dimensionless parameters v and w are related to the 
mass and the elastic constant of the modeli'0, in which the electron is coupled via a harmonic 
force to a fictitious second particle simulating the phonon degrees of freedom. The values of 
V and w can be obtained by the variational approach described by Feynmani and SchultJ. 
In particular, at T = 0, A(g, iS) takes the form: 

_2!« (q^R^ ' 
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The spectral weight follows a Poissonian distribution and it is maximum for an excitation 
involving a number I of phonons of order of / ~ j^. We note that, in all the numerical 
results which will be shown in this paper, the terms with — 6 < / < 6 in the Eq.([TT|) and 
the first six term in the Eq. ([12|) are enough to a good convergence of the real and imaginary 
parts of the correlation function n°(ci;) in the frequency range of interest {u < IOcuq). 

We want to note that the approach proposed in this paper restores, when n —> , the well 
known results of the optical absorption of a single polaronlll and allows to introduce within 
the R.P.A. approximation the effects of the polaron-polaron interaction. This approach is 
expected to be valid when the formation of bipolaron states is not favored [r] = too/ to > -01 
oi 7] < 0.01 and a < 6 — 7)lii. In fact, in this situation the overlap between the wells of 
two particles can be neglected and the system is well approximated by an interacting large 
polaron gas. Of course the proposed approach is valid in the density range where we can 
exclude Wigner-like localization. As discussed by Quemerais et ahlll for coupling constant 
not larger than a = 6 and density n > 10^*^ cm"'^ our approach is justified. 



III. THE RESULTS 

In Fig. 2a is reported the optical absorption per polaron as a function of the frequency for 
different values of the charge carrier density at T = 0. Three different structures appear in 
the normal state conductivity: a) a zero frequency delta function contribution; b) a strong 
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band starting at = cuq that is the overlap of two components: a contribution from the 
intraband process and a peak due to the polaron transition from the ground state to the 
first relaxed excited state; c) a smaller band at higher frequency due to the Frank-Condon 
transition of the polaron. The identification and characterization of these structures stem 
from the detailed analysis of the polaron absorption in the case of non interacting large 
polaronJil. Increasing the charge carrier density, we find that the large polaronic band due 
to the excitation involving the relaxed states ( b contribution) tends to move towards lower 
frequencies while its intensity decreases in favor of the rise of a Drude-like term around 
Lj = (see Fig. 2a and Fig.2b). 

In particular, at T = the Drude weight D, i.e., the coefficient of the zero frequency 
delta function contribution, is determined making use of the sum-rule for the real part of 
the conductivity (see Eq.(^)). This allows an estimate of the effective polaron mass Mf,ff{n) 
as a function of the electron density, being the Drude weight D related to Meff{n) by the 
following expression: 



In table 1 the values of Meff{n) for a = 5 and a = 6 are reported. It is evident that, 
increasing the charge carrier density, the screening of the e-ph interaction increases and the 
polaron mass is reduced, tending to the band mass for large values of n. This behavior is 
a confirmation of the trend obtained at weaker coupling in a model that includes polaron 
screening but neglects exchange effectfll. 

In FigSa we plot the normal state conductivity for larger electron densities at T/T^ = .5, 
where Tq = hujo/kg and kg is the Boltzman constant. Increasing the charge density, we find 
that the optical absorption, in agreement with experimental data in the metallic phase of 
the cupratesEHn^, is more and more controlled by the Drude-like term and, only for very 
high electron densities, no signature of the b) contribution is left . 

The behavior of the polaron absorption with temperature is also of interest. As shown 
in Fig. 3b, with increasing T there is a transfer of spectral weight towards higher frequen- 
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cies. Moreover, the intensity of the contributions b) and c) increases with decreasing the 
temperature saturating at T/Tq ~ 0.5. 

IV. DISCUSSION 

The effects shown are very intriguing in connection of recent measurement of infrared 
absorption in cuprate^l. In fact several experiments on the infrared response of the cuprates 
have pointed out that the optical absorption exhibits features which are common to many 
families of high-T^ superconductors. 

In particular, the infrared normal state conductivity does not diminish with frequency 
as rapidly as one expects from a simple Drude pictureil. This behavior has been interpreted 
in terms of two different models: the anomalous Drude modelii and the multi-component 
or Drude-Lorentz mo del§. In the former approach, the infrared conductivity has been at- 
tributed to a contribution from free charge carriers with an cu-dependent scattering rate 
(arising, for example, by strong interactions with spin wavesii). In the Drude-Lorentz ap- 
proach, instead, one assumes that the absorption is the result of the superposition of different 
structures which can be identified in the conductivity spectra: 1) a Drude-like peak cen- 
tered at zero frequency; 2) infrared-active vibrational modes (IRAV)ii; 3) a broad excitation 
in the infrared band which is constituted by two different components: one temperature- 
independent around .5eV (~ 4OOOcm~^)il0 and the other strongly dependent on temper- 
ature with a peak around .leV (~ 1000 cm~^)i3il ( d band); 4) the charge transfer band 
(CT) in the visible range which is attributed to the charge transfer transitions between 02p 
and Cu3d states. The Drude, the d band and IRAV contributions depend on the charge 
density injected by the doping and a significant transfer of spectral weigth from the d band 
and Irav contributions to the Drude peak is observed increasing the doping. To the d band 
and the IRAV, which appear when extra charge are injected into the lattice of a cuprate, it 
has been assigned a polaronic origin both in electron-doped and hole doped compoundsii'@. 
However, there is no general consensus on the type of polarons involved in the absorption. 
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In particular, the 1000 cm~^ feature has been attributed to optical transitions involving 
smalfl, large polaronsilil or both types of polaronsH. In any case, as first noted by Bi 
and Eklundil in La2-xSrxCu04+s, the polaron tend to be small in the dilute limit (small x) 
whereas in the limit of large x the polaron tend to expand, i.e., the polaron size increases 
with X tracking a dopant-induced transformation to metallic conductivity. 

Finally we note that, recently, Calvani et al.0 have measured the infrared absorption of 
four different perovskite oxides observing an opposing behaviour between the cuprates and 
the non cuprates. In particular, in La2Ni04J^y and in SrMnO^j^y, in which there is a strong 
evidence of the presence of small polarons, the minimum of the d band deepens at low T by 
a transfer of spectral weight towards higher energies. On the contrary, in the slightly and 
heavily doped cuprates at low temperatures the minimum of the d band tends to be filled by 
a transfer of spectral weight in the opposite direction (see also F. Li et al.il,) as it happens 
for a large polaron system previously described. 

The above scenario shows that the normal state conductivity of an interacting large po- 
laron system presents features common to the optical absorption of the cupratesil. Even 
if the similarities found are not conclusive in assigning to large polarons the main role in 
the absorption of cuprates - more accurate calculations are nedeed to attempt a quanti- 
tative comparison between theory and experiments - we believe that the inclusion of long 
range Coulomb interactions among polarons is essential in understanding whether infrared 
absorption in cuprates can be assigned to large polarons. 



V. CONCLUSIONS 

In this paper we calculated the normal state conductivity of a large polaron system in- 
cluding the electrostatic polaron-polaron interaction within the R.P.A. approximation. The 
approach recovers the optical absorption of non interacting Feynman polaron0 and allows 
to introduce in a perturbative way the many polaron effects. With increasing charge carrier 
density, we have found evidence of a transfer of spectral weight of the optical absorption 
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towards lower frequencies: the Drude-like contribution increases as well as the screening 
effects of the e-ph interaction. Moreover, with increasing temperature, there is a transfer of 
spectral weight of the infrared absorption towards higher frequencies. Both these behaviours 
are observed in the infrared spectra of cuprates. 
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TABLE 1 CAPTIONS 

Tab.l. Effective polaron mass, in units of the electron band mass, as a function of the 
charge carrier density. The value of uq is n = 1.4 • 10~^ in terms of R'^, Rp being the 
Frohlich polaron radius. 

FIGURE CAPTIONS 

Fig.l. Diagrammatic representation of the Eq.(|]). The solid line indicates the polaron 
propagator, the phonon is given by the dashed line, the dotted line describes the 
Coulomb interaction and the dotted-dashed line represents the incident photon. 

Fig.2. a) Optical absorption per polaron, at T = 0, as a function of the frequency for 
different electron densities: n = 1.4 • 10~^ (solid line), n = 1.4 • 10~^ (dashed line), 
n = lA ■ 10^^ (dotted line), n = 1.4 ■ 10"^ (dotted-dashed line); b) Optical absorption 
per polaron at finite temperature, T/T^, = .5, for different values of the charge carrier 
density: n = 1.4-10"^ (dashed line), n = 1.4-10"^ (dotted line), n = 1.4 -lO'^ (dotted- 
dashed line). The value of eo/eoo is 3.4. The electron density is measured in units of 
Rp^, Rp being the Frohlich polaron radius, and the conductivity cr{u!) is expressed in 
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terms of ne^ /mujQ. The value of is 7 • 10^° cm ^ when ujq — ?>QmeV and m — nie, 
where rrie is the electron mass. 

Fig. 3. a) Optical absorption per polaron at finite temperature, T/T^ — .5, for different 
values of the charge carrier density: n = 4.5 • 10^^ (sohd hne), n = 1.4 • 10^^ (dashed 
line), n — .28 (dotted line), n = 1.4 (dotted-dashed hne). The electron density is 
measured in units of Rp^, Rp being the Prohlich polaron radius, and the temperature 
is given in units of = huo/kB, kg being the Bohzman constant. In the inset is 
plotted the polaron mobility (/x), i.e. Rea{u! — > 0), as a function of the charge density 
for cuq — SOmeV, eoo = 3 and m — me where m^ is the electron mass; b) Optical 
absorption of a single Feynman polaron for different values of T: T/T^ = (dashed 
line), T/Tu — .5 (dotted line), T /To — 1 (dotted-dashed line). 
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In this paper we calculate the normal state conductivity of a homogeneous three dimensional 
(3D) system. On the other hand it is well known that almost all of the cuprates show strong 
anisotropy and, when the electric field is polarized along the conducting planes, the optical 
response is very close to that of a homogeneous 2D system. In this paper we assume that our 
3D results are qualitatively equal to those of a homogeneous 2D system. This assumption is 
supported by the fact that, for non interacting polarons, a simple scaling relation connects the 
2D and the 3D impedance functions (see F. M. Peeters and J. T. Devreese, Phys. Rev. B 36, 
4442 (1987)). The inclusion of the anisotropy is currently under study and it will be discussed 
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